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Abstract 

Stocking represents the most important management tool worldwide to increase 
and sustain commercial and recreational fisheries in a context of overexploita- 
tion. Genetic impacts of this practice have been investigated in many studies, 
which examined population and individual admixture, but few have investigated 
determinants of these processes. Here, we addressed these questions from the 
geno typing at 19 microsatellite loci of 3341 adult lake trout (Salvelinus namay- 
cush) from 72 unstocked and stocked lakes. Results showed an increase in 
genetic diversity and a twofold decrease in the extent of genetic differentiation 
among stocked populations when compared to unstocked. Stocked populations 
were characterized by significant admixture at both population and individual 
levels. Moreover, levels of admixture in stocked populations were strongly corre- 
lated with stocking intensity and a threshold value of total homogenization 
between source and stocked populations was identified. Our results also suggest 
that under certain scenarios, the genetic impacts of stocking could be of short 
duration. Overall, our study emphasizes the important alteration of the genetic 
integrity of stocked populations and the need to better understand determinants 
of admixture to optimize stocking strategies and to conserve the genetic integrity 
of wild populations. 



Introduction 

Various anthropogenic pressures, such as overexploitation, 
habitats modifications and non-native wild or domesti- 
cated individuals' translocations, may have major impacts 
on the genetic makeup and evolutionary trajectory of wild 
populations, possibly resulting in demographic decline 
(Hutchings 2000; Coltman et al. 2003; Allendorf et al. 
2008; Sharpe and Hendry 2009). This is especially true for 
fish, for which signs of overexploitation through commer- 
cial and recreational fishing have been observed for decades 
(Myers et al. 1997; Hutchings 2001; Post et al. 2002; Cole- 
man et al. 2004; Post 2013). In addition, exploited fish 
populations have also been heavily stocked for decades to 
provide additional harvest opportunities in overexploited 
populations (Waples 1999; Helfman 2007). Although 
stocking is a very important tool to meet management 
goals, it may negatively impact the evolutionary potential 



of populations (reviewed in Araki et al. 2008; Fraser 2008; 
McClure et al. 2008) as well as their genetic integrity 
(reviewed in Laikre et al. 2010). 

Mounting evidence indicates that captive reared fish 
released in the wild can experience a lower fitness than 
their wild conspecifics, even after very few generations of 
captivity (Ford 2002; McGinnity et al. 2003; Araki et al. 
2007a, 2008; Theriault et al. 2011; Christie et al. 2012; 
Milot et al. 2013; but see Berejikian et al. 2009). Domesti- 
cation selection seems the most likely explanation for this 
fitness reduction (reviewed in Araki et al. 2008) as captive 
reared fish experience different selection regimes than wild 
fish (Heath et al. 2003; Blanchet et al. 2008; Fraser 2008; 
Williams and Hoffman 2009; Christie et al. 2012). In addi- 
tion, captive reared fish potentially adapted to captivity 
and distinct from the wild populations can lower stocked 
population fitness through introgressive hybridization 
(McGinnity et al. 2003; Araki et al. 2008, 2009). Outbreed- 
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ing depression may then occur if the interbred offspring 
suffers a disruption of coadapted genes complex and/or a 
disruption of the interactions between genes and the envi- 
ronment they evolved in (reviewed in Allendorf et al. 
2001; reviewed in Edmands 2007; Bougas et al. 2010; 
Granier et al. 2011). 

Modification of neutral patterns of genetic diversity is 
also often associated with stocking. Stocked populations 
may show a loss of genetic diversity (Eldridge et al. 2009), 
a loss of genetic differentiation with other populations 
(Susnik et al. 2004; Eldridge and Naish 2007; Eldridge et al. 
2009; Hansen et al. 2010; Marie et al. 2010; Lamaze et al. 
2012; Perrier et al. 2013b) and a displacement of the local 
gene pool (reviewed in Laikre et al. 2010). Genetic integrity 
of locally adapted populations can also be threatened 
through introgressive hybridization (McGinnity et al. 
2003, 2009; Bourret et al. 2011). Numerous studies docu- 
mented various degrees of population and individual 
admixture in wild populations subjected to stocking (Sus- 
nik et al. 2004; Hansen et al. 2009; Marie et al. 2010; 
Dawnay et al. 2011; Karaiskou et al. 2011; Perrier et al. 
2011; Lamaze et al. 2012). However, the link between 
stocking intensity and the extent of modification of local 
genetic diversity, including allelic richness and admixture, 
have been much less investigated (Hansen 2002; Marie 
et al. 2010; Perrier et al. 2011, 2013b; Lamaze et al. 2012). 
Moreover, despite the importance of such knowledge for 
population restoration, few studies attempted to predict 
the evolution of the genetic composition of stocked popu- 
lations when stocking has ceased (but see Hansen and 
Mensberg 2009; Hansen et al. 2010; Perrier et al. 2013b). 

Here, we assess the impact of stocking intensity on the 
evolution of the genetic composition of a large number of 
lake trout (Salvelinus namaycush, Walbaum 1792) popula- 
tions from Quebec, Canada. Lake trout is a freshwater fish 
native to North America and inhabiting cold-water lakes. It 
is long living (up to 49 years) and becomes sexually mature 
between 7 and 13 years and is iteroparous (Scott and 
Crossman 1998). A pronounced genetic structure has been 
documented among wild populations (Piller et al. 2005; 
Halbisen and Wilson 2009; Northrup et al. 2010; McCrac- 
ken et al. 2013). Also, Piller et al. (2005) and Halbisen and 
Wilson (2009) reported noticeable impacts of stocking on 
the genetic integrity of some stocked populations. In the 
province of Quebec, Canada, numerous lakes were stocked 
for over 40 years (about five generations) and the stocking 
history was exhaustively recorded, thus providing a unique 
context for documenting genetic changes of wild popula- 
tions as a function of stocking history. Here, we specifically 
aim to (i) document genetic diversity and structure among 
unstocked lake trout populations, (ii) quantify the impacts 
of stocking on the distribution of genetic diversity within 
and among stocked lake trout populations, (iii) determine 
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the influence of stocking intensity on admixture rates 
within stocked populations and on homogenization among 
source and targeted populations, and to (iv) investigate the 
evolution of the genetic composition of stocked popula- 
tions after stocking has ceased. 

Materials and methods 

Study site, sample collection and stocking 

A total of 3341 adult lake trout were sampled in 72 lakes 
from 10 independently managed and stocked administra- 
tive regions across the province of Quebec, Canada (Fig. 1; 
Table 1, Table S4). An exception was the Chaudiere-Appal- 
aches and Estrie regions, which were both stocked using 
the same broodstocks. Therefore, lakes from these two 
regions were analyzed as if they belonged to the same 
region (Estrie). Adipose or pelvic fin clips were collected 
between 1997 and 2011 for 3003 of those individuals with 
the help of anglers and the provincial government. These 
samples were preserved in 95% ethanol. Dry scales of 338 
adult lake trout from four out of 72 lakes were also col- 
lected between 1974 and 1983 (Table 1). The mean number 
of individuals sampled per lake was 42 (range 15-68). 

The stocking history of each stocked lake is recorded by 
the provincial government since 1900 in a central database. 
For each stocking event, a standardized stocking form is 
filed, compiled in the database and carefully verified to 
detect potential errors. Although the database is reliable, 
we performed additional validation of the data for the 
stocked lakes used in this study for further improvements. 
Lakes were categorized as 'unstocked' when no stocking 
was recorded and as 'stocked' when at least one stocking 
event was documented. The mean number of stocking 
events in stocked lakes was 18 (range 5-50). A source lake 
was a lake where broodstock was collected and used to pro- 
duce juveniles in captivity, which were released in wild 
populations, usually at age 1+. Of the 72 sampled lakes, 10 
were source lakes (n = 432), 40 were unstocked lakes 
(n = 1600) and 22 were stocked lakes (n = 971). Of the 10 
source lakes, two were never stocked, whereas eight were 
stocked with local and/or nonlocal individuals. Among 
lakes sampled between 1974 and 1983, three were never 
stocked and one was stocked with juveniles produced with 
its own broodstock after 1962. These samples were used to 
validate the temporal stability of genetic characteristics 
within populations never stocked with non-native individ- 
uals as contemporary samples were also collected as 
described above. 

DNA extraction 

DNA was extracted from adipose fin and pelvic fin clip 
(2 mm 2 ) using a modified version of Aljanabi and 
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Figure 1 Geographical locations of sampled populations in the province of Quebec, Canada (see Table 1 for more information). 



Martinez (1997) salt extraction method. DNA precipita- 
tion was performed with isopropanol for 30 min. After 
washing the pellets, we performed a centrifugation at 
10 000 g for 10 min at 4°C. Finally, we eluted the sam- 



ples in 100 fi~L sterile dH 2 0. DNA from dried scales was 
extracted using the Qiagen DNeasy Blood and Tissue kit 
(Qiagen inc., Valencia, CA, USA) following a modified 
version of the manufacturer's Bench Protocol for animal 
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tissues (spin-column protocol). Namely, we incubated 
samples with proteinase K at 37°C for 3 days. After hav- 
ing discarded flow-through (buffer AW2), we performed 
a second centrifugation for 1 min at 20 000 g. We then 
heated up buffer AE at 56°C and used 75 /(L of it to 
perform a first elution. We repeated the elution with a 
second 75 fiL warm AE buffer. 

DNA amplification and microsatellite genotyping 

Individuals collected between 1997 and 2011 (n = 3003) 
were analyzed at 19 microsatellite loci using Qiagen Multi- 
plex PCR Kit (Qiagen inc.) (Table SI, Supporting informa- 
tion). Each 10 fiL multiplex PCR reaction contained 
5-25 ng DNA template, 5 /iL Qiagen multiplex reaction 
buffer and forward and reverse primers at different concen- 
trations: multiplex A contained 4.5 fiM SfoD75, 3 ftu 
Sfo308lav, 3.5 fm Sfo226lav, 0,8 fm SnaMSU02, 1 fm 
SnaMSU03; multiplex B contained 3 fm SnaMSU06, 
0.2 fm SnaMSU08, 6.5 fm SnaMSU09, 7.0 fm SnaMSUlO, 
7 fm Sco202; multiplex C contained 0.5 fm SnaMSUll, 
2 fiu SnaMSUU, 1.5 fm SnaMSUH, 10 fiM SnaMSU07, 
0.25 fm Sco215; multiplex D contained 1 fiM SnaMSUOl, 
2.5 fiM Sco200, 0.3 fiM Smm22, 1.2 fm Sssp2201 (Crane 
et al. 2004; Dehaan and Ardren 2005; Paterson et al. 2004; 
Perry et al. 2005; Rollins et al. 2009; T. L. King, S. E. Julian, 
R. L. Coleman and M. K. Burnham-Curtis, unpublished). 
Amplifications were performed using a Tl Biometra ther- 
mocycler (Biometra, Kirkland, QC, Canada) with a 15 min 
activation step at 95°C, followed by 35 cycles of denatur- 
ation at 95°C for 30 s, annealing for 3 min at 60°C (multi- 
plex A, B and C) or 56°C (multiplex D), extension at 72°C 
for 1 min, followed by a final extension at 72°C for 
10 min. Individuals collected between 1974 and 1983 were 
genotyped with the same loci except for Sfo226Lav, Sfo308- 
Lav, Sco215, SnaMSU07 and Sssp2201 which did not 
amplify or amplified poorly with old DNA template. These 
individuals were amplified in multiplex with the same PCR 
programs except for multiplex A for which we used an 
annealing temperature of 56°C instead of 60°C. Amplified 
products were migrated via electrophoresis using an ABI 
3130xZ capillary DNA sequencer (Applied Biosystems Inc.). 
Genescan 500 LIZ (Applied Biosystems Inc., Burlington, 
ON, Canada) was used as a standard to determine allele 
sizes, which were scored using genemapper 4.0 (Applied 
Biosystems Inc). 

Genetic diversity 

The software micro-checker (Van Oosterhout et al. 2004) 
was used to assess the potential presence of null alleles and 
large allelic dropout. P IS and number of alleles per locus 
(A) were estimated using fstat 2.9.3 (Goudet 2001). Allelic 
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richness (A R ) adjusted for the smallest sample size (n = 15) 
was calculated using HP-Rare 1.0 (Kalinowski 2005). We 
calculated mean A R for unstocked and stocked populations 
and tested for difference using Student's f-test. We used a 
linear model to assess the relationship between A R and the 
natural logarithm of lake area [ln(LakeArea)] for unstocked 
and stocked lakes. Strength of the correlation and signifi- 
cance were assessed with Pearson's product-moment corre- 
lation coefficient. Expected and observed heterozygosities 
(H E and H Q respectively) were obtained with genetix 4.05 
(Belkhir et al. 1996-2004). We calculated mean H E for 
unstocked and stocked populations and tested for differ- 
ence using Student's r-test. Effective population size (N E ) 
was estimated using the linkage disequilibrium method 
implemented in LDNe (Waples and Do 2010) for 
unstocked lakes only. N E would have been meaningless in 
stocked lakes due to biases resulting from individual and 
population admixture which increase linkage disequilib- 
rium (Araki et al. 2007b). We used a linear model to assess 
the relationship between N E and ln(LakeArea) and between 
N E and A R . Strength of the correlations and significance 
were assessed with Pearson's product-moment correlation 
coefficient. 

Population genetic structure 

We quantified the extent of genetic differentiation between 
each pair of populations using pairwise F ST measures with 
fstat 2.9.3 (Goudet 2001). We tested for difference in 
genetic differentiation between the unstocked versus 
stocked lakes with fstat using the comparison among 
groups with 10 000 permutations. To assess the relation- 
ship between stocking intensity and genetic differentiation 
among source and stocked populations, we used Spear- 
man's rank test correlation. Stocking intensity was calcu- 
lated as the total number of fish stocked per hectare and as 
the number of stocking events (years). 

We assessed the distribution of genetic diversity among 
administrative regions, among populations within regions, 
and within populations using analyses of molecular vari- 
ance (amova) implemented in arlequin 3.5 (Excoffier and 
Lischer 2010) and tested with 10 000 permutations. To 
assess the influence of stocking on the distribution of 
genetic diversity among and within regions, three different 
amovas were performed. The first amova was performed 
with unstocked populations (n = 42) grouped by 10 
regions (Abitibi-Temiscamingue, Bas-Saint-Laurent, 
Chaudiere-Appalaches-Estrie, Cote-Nord, Lanaudiere, Lau- 
rentides, Mauricie, Nord-du-Quebec, Outaouais, Sague- 
nay- Lac-Saint-Jean). The second amova was performed 
with stocked populations (n = 30) grouped by a subset of 
five regions where stocking occurred (Bas-Saint-Laurent, 
Chaudiere-Appalaches-Estrie, Lanaudiere, Laurentides, 
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Outaouais) as the other five did not contain any stocked 
populations. The third amova was performed with these 
same five regions but with only the unstocked populations 
[n = 16) for strict comparison with the above-mentioned 

AMOVA. 

Admixture analyses 

Individual clustering was achieved using the Bayesian 
method implemented in structure 2.3.1 software (Prit- 
chard et al. 2000). This analysis was carried out to assess 
the number of genetic clusters present in our data and 
quantify admixture proportions at both the individual and 
population levels among stocked and source populations. 
structure analyses were performed assuming an admixture 
model without priors. We ran the analysis for seven differ- 
ent lake groupings (defined below) with k genetic clusters 
from 1 to the number of lakes plus 3 (number of lakes var- 
ied from 9 to 42 according to the group of lakes consid- 
ered), and with 15 replicates for each k. The first group was 
composed of all unstocked lakes only. The second to the 
sixth groups each represented one of the five stocked 
administrative regions, composed of unstocked lakes, 
stocked lakes, and all known source lakes used for stocking 
within each region. The last group was composed of the 
four lakes with historical and contemporary samples. 

Each Structure run started with a burn-in period of 
50 000 followed by 300 000 Markov Chain Monte Carlo 
(MCMC) steps. For each lake group considered, we 
selected the best k according to the variation of likeli- 
hood (Pritchard et al. 2000), the AA; method (Evanno 
et al. 2005) implemented in structure harvester (Dent 
and vonHoldt 2012) and consistency of structure out- 
puts. The software distruct (Rosenberg 2004) was used 
to plot structure outputs. 

To further investigate the impact of stocking, we quanti- 
fied proportions of fish belonging to local clusters, nonlocal 
clusters, or to a category of individuals putatively admixed 
between local and nonlocal populations and representing 
putative hybrids. For this analysis, we used lakes sampled 
within the five stocked regions and excluded source lakes. 
Following Vaha and Primmer (2006), individuals were clas- 
sified as local, nonlocal and putatively admixed when their 
individual admixture proportion (ind. q-values) of the local 
cluster were, respectively, >0.70, <0.30, and [0.30 and 
0.70]. For statistical analyses, stocked and unstocked lakes 
were separated into distinct groups. A chi-square test was 
then used to test whether differences in proportions of 
local, admixed, and nonlocal individuals between the 
unstocked and the stocked groups were significant. This 
test was repeated for unstocked and stocked groups within 
each of the five stocked regions. Chi-square was also used 
to test whether differences in the proportion of these indi- 



viduals were significant among stocked as well as unstocked 
groups between regions. 

Effects of stocking intensity on introgressive hybridization 

We investigated the effects of stocking intensity on popula- 
tion admixture (q-membership) in two ways using R [R 
Development Core Team (2012)]. Firstly, we applied a 
logistic regression model to investigate the link between 
stocking intensity and the membership of stocked popula- 
tions to the local cluster (local q-membership). Unstocked 
populations were included in the regression as they were 
used to assess local q-membership when no stocking was 
documented. We attributed a value of '1' to populations 
with a local q-membership >0.50 and a value of '0' to popu- 
lations with a local ^-membership <0.50 as few intermedi- 
ate values were found (see Results section). Secondly, we 
examined the link between the number of stocked fish per 
hectare and the membership of stocked populations to 
their source(s) cluster {source q-membership) as well as the 
link between the number of stocking events and the source 
q-membership using Spearman's rank test correlations. For 
these analyses, source q-membership was used, as many 
sources could have been used for each stocked lake. 

Lastly, we assessed the effect of the time spent (in years) 
since the last stocking event on the populations' source 
q-membership. For this latter analysis, as stocking pressure 
may influences populations' source q-membership, we only 
used pairs of source and stocked populations that 
exchanged similar number of fish per hectare to eliminate 
potential bias originating from differences in stocking 
intensity. The only stocking intensity in number of fish 
stocked per hectare for which we had consistent data over 
the 40 years period was 17-22 fish stocked per hectare. 
Nine pairs of source and stocked populations corresponded 
to this criterion (# unstocked- # source: 11-47, 11-68, 
14-10, 20-68, 22-21, 23-68, 32-48, 33-68, 42-37). For this 
analysis, we could not find pairs of source and stocked 
populations with similar number of stocking events over 
the 40 years period. We used Spearman's rank test correla- 
tion to assess significance of the correlations. 

Results 

Genetic diversity 

Population genetic diversity indices for each population are 
shown in Table 1. Among the contemporary samples (col- 
lected between 1997 and 2011), 72 of 1368 permutation 
tests conducted in fstat revealed significant F IS values. 
However, micro-checker showed that only seven of these 
72 significant P IS were associated with null alleles. At the 
population level, six of the 72 P IS tests were significant. Of 
these, one was associated with the potential presence of null 
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alleles (#71) and two corresponded to stocked populations 
(#10 and 11). For the historical samples (collected between 
1974 and 1983), six of 70 permutation tests yielded signifi- 
cant F 1S values of which five were associated with null 
alleles for five different loci. At the population level, one of 
five F IS was significant and associated with the potential 
presence of null alleles (#59). 

All 19 loci were moderately to highly polymorphic 
with the number of alleles per locus ranging from 8 to 
56 with an average of 32. Over all loci and populations, 
a total of 613 alleles were identified. Global A R was 
11.96 over all loci and populations, while mean A R per 
population was 7.29 and ranged from 3.00 to 11.54. A R 
was significantly higher in stocked (mean A R = 8.86 
ranging from 5.77 to 11.15) compared unstocked popu- 
lations (mean A R = 6.16 ranging from 3.00 to 11.54) 
(t = -6.15; d.f. = 57.61; P < 0.0001; Fig 2A). The In 
(LakeArea) was positively correlated with A R among 
unstocked populations (R 2 = 0.41; t= 5.44; P < 0.0001), 
but not among stocked populations (R 2 = 0.03; 
t = 0.257; P = 0.80). Average H E was 0.695 over all pop- 
ulations ranging from 0.415 to 0.856. As observed for 
A R , H E was significantly higher among stocked popula- 
tions (mean H E = 0.774 ranging from 0.651 to 0.856) 
than among unstocked ones (mean H E = 0.639 ranging 
from 0.4147 to 0.836) (f = -6.71; d.f. = 64.11; 
P < 0.0001). A/g varied from 2 to 2509 among unstocked 
populations (with a median value of 134). We detected 
a significant positive correlation between N E and lake In 



(LakeArea) (_R 2 = 0.26; t = 3.92; P < 0.001), and between 
N E and A R (R 2 = 0.09; f = 2.25; P = 0.0298). 

Population genetic structure 

Results showed that stocking had a profound effect on pat- 
terns of population structure. Estimates of global _F S t 
among the 42 unstocked populations was more than twice 
of that observed among the 30 stocked populations {0.227, 
(CI [0.197-0.260]) vs. 0.102 (CI [0.082-0.123]), 
P < 0.0001}. Pairwise F ST values ranged from 0.025 to 
0.531 among unstocked populations and from —0.016 to 
0.298 among stocked populations (Fig 2B; Table S2, Sup- 
plementary material). Among unstocked populations, the 
variance of pairwise -Fst values was large with a median 
value of 0.246, but much smaller among stocked popula- 
tions with a median value of 0.088 (Fig. 2B). Moreover, a 
highly significant negative correlation was observed 
between pairwise F ST values comparing source versus 
stocked populations and stocking intensity, either in num- 
ber of fish stocked per hectare (p = -0.432; P < 0.0001; 
Fig 3A) or in number of stocking events (p = —0.489; 
P < 0.0001; Fig 3B). amovas also revealed that genetic dif- 
ferentiation among populations within region was much 
more pronounced for unstocked populations with 
P sc = 0.199 (CI [0.173-0.226]) than stocked F sc = 0.041 
(CI [0.033-0.050]) (P< 0.0001; Table 2). In contrast, 
genetic differentiation among regions was twice lower for 
unstocked populations with F C t = 0.036 (CI [0.026- 
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Figure 2 (A) Distribution of A R value frequency for unstocked (N populations = 42; white bars) and stocked (W populations = 30; dark gray bars) 
populations; triangles indicate median A R value for each distribution; light gray is the overlapping between the two distributions, (B) distribution of 
pairwise F S t frequency among unstocked (N pairwise comparisons = 861; white bars) and among stocked populations (N pairwise compari- 
sons = 435; dark gray bars); triangles indicate median F S t for each distribution; light gray is the overlapping between the two distributions. 
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Figure 3 (A) Pairwise F S t values between source and stocked populations (N pairwise comparisons = 78) as a function of stocking intensity in num- 
ber of fish per hectare and B) in number of stocking event with a 0.95 confidence envelope, C) populations' membership (N pop. = 38; N 
ind. = 1629) to their own genetic cluster (local q-membership) as a function of stocking intensity in number of fish stocked per hectare and D) in 
number of stocking events, E) stocked populations' membership to their source populations cluster (source q-membership) as a function of number 
of fish stocked per hectare (colors) (N pairwise comparisons = 69), with a 0.95 confidence envelope. 
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Table 2. Global genetic differentiation between populations (F ST ) as calculated by fstat, among lakes within regions structure (F sc ) and among 
regions structure (F CT ; 42 populations within 1 0 regions for the all unstocked populations group, 1 6 populations within five regions for the subset of 
unstocked populations group and 30 populations within five regions for the stocked group) as calculated by Arlequin. 





F ST [5-95%] 


F sc [5-95%] 


F ct [5-95%] 


All unstocked populations 


0.227 [0.197-0.260] 


0.199 [0.173-0.226] 


0.036 [0.026-0.048] 


Subset unstocked populations 


0.184 [0.155-0.221] 


0.166 [0.143-0.191] 


0.029 [0.017-0.042] 


Stocked populations 


0.102 [0.082-0.123] 


0.041 [0.033-0.050] 


0.064 [0.050-0.078] 



0.048]) as compared to F CT = 0.064 (CI [0.050-0.078] for 
stocked (P < 0.0001; Table 2). 

Pariwise F ST values among historical and contemporary 
populations were not different from 0 except for popula- 
tion #59 for which a small F ST value of 0.01 was found 
(Table S3). 

Admixture analyses 

Most genetic clusters delineated by structure corresponded 
to a unique unstocked lake. Among the 42 unstocked pop- 
ulations, k = 39 was the most likely number of clusters 
(Fig 4A). In six cases, one cluster corresponded to two 
lakes. Also, two lakes comprised more than one cluster. 
However, performing independent clustering analyses for 
each region, each unstocked population corresponded to a 
single cluster except lakes #25 and #26, which corresponded 
to a same cluster, and lakes #44 and #47, which corre- 
sponded respectively to two and four clusters. Lakes #25 
and #26 were in fact geographically close to each other 
(2.7 km) and physically connected. Considering both 
unstocked and stocked populations on a regional scale, the 
best clustering solutions were k = 8 for the Bas-Saint-Lau- 
rent region, k = 6 for the Estrie-Chaudiere-Appalaches 
region, k = 6 for the Lanaudiere region, k = 17 for the Lau- 
rentides region and k = 1 1 for the Outaouais region 
(Fig 4B). In every stocked region, stocked populations 
shared various proportions of ancestry with their source(s), 
while unstocked populations corresponded to a unique 
cluster (except lakes 25 and 26; Fig. 4B; Table 3). Finally, 
k = 4 was retained in the case of the comparison of histori- 
cal and contemporary samples (Fig 4C). Each unstocked 
population corresponded to the same unique cluster in 
both contemporary and historical samples. 

Admixtures analyses also revealed that stocking had a 
profound impact on the genetic integrity of lake trout pop- 
ulations. Considering the whole dataset, the proportion of 
fish assigned to the local cluster, a nonlocal cluster or a 
putatively admixed category differed between unstocked 
and stocked lakes (x 2 {2) = 991, P < 0.0001; Fig 5A). Over 
the 16 unstocked populations sampled in the five different 
stocked regions, 97.01% ± 0.99% of the fish belonged to 
local cluster on average while no fish belonged to nonlocal 
cluster and 2.99% ± 0.99% of the fish belonged to the 



putatively admixed category (Fig. 5A). In contrast, among 
the 22 stocked populations, 17.49% ± 6.61% of the fish 
belonged to the local cluster, 56.48% ± 7.06% of the fish 
belonged to a nonlocal cluster, and 26.02% ± 3.78% of the 
fish were putatively admixed (Fig. 5A). Such differences 
were also observed within each region and were highly sig- 
nificant (Bas-Saint-Laurent: x 2 (i) = 258, P < 0.0001; Estrie: 
X 2 ( 2 ) = 163, P < 0.0001; Lanaudiere: y\ 2) = 159, 
P < 0.0001; Laurentides: x 2 {2) = 245, P < 0.0001 and Ou- 
taouais: x 2 (2) = 213, P < 0.0001; Fig 5B-F). While we did 
not observe any significant difference in the proportion of 
fish assigned to the local cluster, nonlocal cluster or puta- 
tively admixed category among the five regions for 
unstocked populations (x 2 (&) = 8.66, P = 0.372), we did 
observe differences among regions for the stocked popula- 
tions (x 2 (s) = 311, P < 0.0001). In particular, while almost 
all individuals were assigned to a nonlocal cluster in 
stocked populations from Estrie, Lanaudiere and Bas- 
Saint-Laurent regions, a noticeable proportion of local fish 
were detected in stocked populations from Outaouais and 
Laurentides regions. 

Effects of stocking intensity 

Logistic regressions showed that populations' membership 
to local cluster sharply decreased with the increase of stock- 
ing intensity (Fig. 3C,D). With a small stocking intensity of 
eight fish per hectare (Fig. 3C) or after a single stocking 
event (Fig. 3D), the logistic regression predicted that 10% 
of the stocked populations show a 0.00 local q-membership 
value to the local cluster. For a stocking intensity of more 
than 86 fish per hectare or more than 18 stocking events, 
the regression predicted that 100% of the stocked popula- 
tions show a 0.00 local q-membership value to the local 
cluster. For each pair of source and stocked population, a 
strong positive correlation was observed between the stock- 
ing intensity and the proportion of membership of stocked 
populations corresponding to their source population 
(source q-membership) (Fig. 3E; p = 0.697; P < 0.0001). 
With <20 fish stocked per hectare, the contributions of 
source populations were small (source q-member- 
ship = 0.05 in average). According to the locally weighted 
scatterplot smoothing (LOESS), the source q-membership 
reached 0.5 with an average of 45 fish per hectare being 
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Figure 4 Bayesian individual clustering results with structure for A) unstocked populations (k = 39); (B) regions where stocking occurred, with 
the unstocked populations on the left, the stocked populations in the middle and the source populations used for each region on the right 
(Bas-Saint-Laurent for k = 8, Estrie-Chaudiere-Appalaches for k = 6, Lanaudiere for k = 6, Laurentides for k = 14 and Outaouais for k= 11); 
(C) historical samples compared to contemporary samples for k = 4. Colored columns represent proportions of membership of each individual 
to each cluster. Green circles represent unstocked populations, downward pointing red triangles represent stocked populations and blue trian- 
gles represent source populations. 
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Table 3. Membership of stocked populations attributed to their stock- 
ing source lakes (source g-membership). Underlined values represent 
putative undocumented stocking between lakes while italic data repre- 
sent putative indirect admixture via other source lakes that were them- 
selves stocked. Dashes indicate that stocking was not documented in 
our analyses. 

Source lake no 

Stocked 

lake no 7 10 21 36 37 40 47 49 68 

11 - 0.61 0.10 - - 0.03 

12 - 0.75 0.10 - - 0.05 

13 - 0.56 0.11 - - 0.25 

14 - 0.87 0.06 - - 0.01 
20 0.01 - 0.88 - - 0.01 0.03 

22 0.01 0.93 - 0.03 0.02 

23 0.00 - 0.97 - - 0.01 0.01 

29 0.03 0.45 0.03 0.02 

30 0.10 O04 0.60 - 0.10 

31 0.01 - 0.07 0.08 0.14 - 0.01 0.01 0.11 

32 - 0_]6 - - 0.01 0.04 

33 0.15 042 0.45 - 0.09 0.01 0.09 

34 O03 - - 0.02 

35 0.21 O03 0.35 - 0.05 0.03 0.21 

41 0.01 0.95 - 0.01 0.04 

42 - 0.24 0.68 - - 0.05 
56 - 0.91 

63 - - 0.51 

64 - 0.09 

65 - 0.87 

66 - 0.87 

67 ________ 0.10 



stocked. However, the contribution of source populations 
was very variable depending on the pairs of stocked and 
source populations. 

Finally, there was a negative correlation between time 
elapsed since last stocking event and stocked populations' 
source q-membership, for each pair of source and stocked 
population (p = -0.854; P < 0.01; Fig. 6). The LOESS 
regression predicted that the source q-membership 
decreased from more than 0.95 to 0.20 approximately 
15 years following the last stocking event, and that the con- 
tribution of the source population was almost null (source 
q-membership < 0.10) 18 years following the last stocking 
event. 



Discussion 

Our main objective was to investigate the impacts of stock- 
ing on the distribution of genetic diversity within and 
among wild lake trout populations over a very large geo- 
graphic range. Our results demonstrate that stocking 
resulted in an increase in neutral genetic diversity in 
stocked populations compared with unstocked ones. Given 
that stocking practices has traditionally been executed at 



the administrative regional scale, we observed a decrease in 
differentiation among populations within regions and an 
increase in regional genetic differences due to stocking. 
While unstocked populations were largely composed of 
local and nonadmixed individuals, stocked populations 
exhibited variable admixture rates with stocking sources. 
Furthermore, this study provides a rare clear link between 
stocking intensity and admixture rates (but see Marie et al. 
2010). Finally, our results suggest a decrease over time of 
admixture in stocked populations after stocking has ceased. 

Effect of stocking on allelic richness in lake trout 
populations 

In line with several previous studies (Ihssen et al. 1988; 
Halbisen and Wilson 2009; Northrup et al. 2010), genetic 
diversity (A R ) in unstocked lakes was strongly correlated 
with lake size (R 2 = 0.41). As lake trout carrying capacity is 
proportional to the lake area (Shuter et al. 1998), as is the 
effective population size (N E ; this study and McCracken 
et al. 2013), variation in population size and N E among 
lakes probably influences the extent of population genetic 
diversity through variable intensity of genetic drift 
(Charlesworth 2009; Frankham et al. 2010). Interestingly, 
we found that the average local genetic diversity 
(A R = 6.16) was more than two times lower than the global 
diversity (A R = 11.75), which is typical of a nonmigratory 
species where populations are isolated from one another 
(Waples and Gaggiotti 2006). Overall, genetic drift may 
have been a determinant process in the erosion of the origi- 
nal genetic diversity, especially given the probable absence 
or very limited potential for contemporary dispersal among 
lakes since the end of colonization (6000 years ago; Wilson 
and Hebert 1996). 

The increase of genetic diversity in stocked compared 
with unstocked lakes as well as the breakdown of the corre- 
lation between genetic diversity and lake area for stocked 
populations reflects the introduction of new alleles from 
nonlocal populations through stocking. Similar results were 
reported for brook char (S. fontinalis) where stocked popu- 
lations with a domestic strain also showed a higher genetic 
diversity than unstocked ones (Marie et al. 2010). In that 
particular case, stocked populations shared private alleles 
with nonlocal domestic broodstock used for stocking. 
Interestingly, it has been shown that when wild local 
broodstock were used as sources of stocking, no or little 
difference was observed in the level of genetic diversity fol- 
lowing stocking (Eldridge and Killebrew 2008; Small et al. 
2009; Gow et al. 2011). In our study, the local increase in 
genetic diversity was not associated with an increase in the 
global genetic diversity (A R = 11.75 among all unstocked 
populations, A R = 11.05 among the subset of unstocked 
populations from stocked regions only, A R = 11.35 among 



636 



©201 4 The Authors. Evolutionary Applications published by John Wiley & Sons Ltd 7 (2014) 625-644 



Valiquette et al. 



Genetic impacts of stocking in lake trout 



(A) All populations (B) Bas-Saint-Laurent (C) Estrie 




Unstocked Stocked Unstocked Stocked Unstocked Stocked 

Figure 5 Percentages (±SE) of fish assigned to local cluster, putatively admixed category (put. Admixed) or nonlocal cluster for unstocked and 
stocked populations: (A) all populations sampled in stocked regions (N pop. = 38; N ind. = 1629), B) Bas-Saint-Laurent (N pop. = 6; N ind. = 258), 
C) Estrie-Chaudiere-Appalaches (N pop. = 4; N ind. = 174), D) Outaouais (N pop. = 1 1; N ind. = 504), E) Lanaudiere (N pop. = 4; N ind. = 159), F) 
Laurentides (N pop. = 12; N ind. = 534). 



genetic diversity of wild stocked populations by increasing 
local allelic richness through the addition of nonlocal 
alleles. 

Effect of stocking on genetic structure among populations 

Overall, the genetic differentiation among unstocked popu- 
lations was pronounced (global f ST = 0.23) and consistent 
with other studies on lake trout (Halbisen and Wilson 
2009; Northrup et al. 2010; McCracken et al. 2013) as well 
as with values observed in other nonmigratory salmonids 
(global F sr brook charr 0.23, Marie et al. 2010; bull trout 
0.33, Taylor et al. 2001; freshwater Atlantic salmon 0.14, 
Tessier et al. 1997; European grayling 0.24, Koskinen et al. 
2002). Such large genetic differentiation is most probably 
linked to significant genetic drift and absence or very lim- 
ited gene flow among populations, as discussed above. In 
contrast, the twofold decrease in genetic differentiation 
found among stocked populations (global F sr 0.10) com- 
pared with unstocked populations indicates a strong 
genetic homogenization effect induced by stocking at the 
regional scale. This tendency of genetic homogenization 
among stocked populations has also been observed in 




Number of years since last stocking event 

Figure 6 Stocked populations' membership to their source populations 
cluster (source g-membership) as a function of time elapsed (in years) 
since the last stocking event (N pop. = 9), with a 0.95 confidence enve- 
lope. 

stocked populations), further indicating a homogenization 
of the genetic diversity among stocked populations rather 
than an increase of the global genetic diversity. These 
results overall confirm that stocking fish from wild but 
nonlocal origins can significantly modify the extent of 
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previous studies in several salmonids (e.g., brook charr: 
Marie et al. 2010; Lamaze et al. 2012; lake trout: Halbisen 
and Wilson 2009) and is due in all cases to the broad scale 
practice of using few source populations for broodstock. 

Our study also revealed that the genetic differentiation 
between source and stocked populations decreased with 
stocking intensity. Reduction of the genetic differentiation 
between source and stocked population was documented 
in European grayling Thymallus thymallus (Koskinen et al. 
2002) and Atlantic salmon Salmo salar (Perrier et al. 2011, 
2013b) and has been linked to the stocking history and the 
presence of nonlocal fish and/or putatively admixed indi- 
viduals. However, these studies did not empirically investi- 
gate how stocking intensity may have influenced levels of 
genetic differentiation between source and stocked popula- 
tions. Here, we observed that while relatively low-stocking 
intensities ( 1-20 fish/ha or 1-4 stocking events) resulted in 
highly variable levels of genetic differentiation between 
source and stocked populations, high stocking intensities 
(more than 40 fish/ha or more than eight stocking events) 
resulted in a nearly zero genetic differentiation between 
these groups. To our knowledge, this is the first attempt to 
empirically determine a threshold value of total homogeni- 
zation between source and stocked populations. Such a 
threshold obviously has several implications in terms of 
management of wild populations as discussed below. 

Artificial hierarchical structuring induced by regional 
stocking 

The weak hierarchical structuring observed among 
unstocked populations most probably reflects the already 
evocated influence of genetic drift and the absence of 
migration among populations leading to a large and not 
geographically hierarchical differentiation of populations. 
This observation is concordant with large-scale phylogeog- 
raphy studies of lake trout mtDNA which suggest little con- 
temporary dispersion among populations from distinct 
mitochondrial lineages since the end of colonization 
6000 years ago (Wilson and Hebert 1996, 1998). Besides, 
Quebec's territory was mainly colonized by a single lineage 
from the Atlantic refuge (Wilson and Hebert 1996) which 
may have limited potential lineages that could have 
resulted in hierarchical structuring. Other studies con- 
ducted on lake trout also documented little or no evidence 
of geographical grouping (Halbisen 2008; McCracken et al. 
2013). Thus, each unstocked lake trout population could 
be considered as a unique 'genetic entity' without any clear 
regional grouping. 

Considering only stocked populations, a more pro- 
nounced genetic differentiation was observed among 
administrative regions than among lakes within regions, in 
total opposition with results obtained with unstocked lakes. 



Valiquette et al. 

More precisely, while the genetic variation among stocked 
populations was four times lower than unstocked, it was 
two times higher among regions. These results are coherent 
with the managing strategies used for the last 20 years 
where stocking has been conducted using one or two 
broodstock sources to stock lakes from each region. Such a 
regionally based stocking practice promoted the homogeni- 
zation of the genetic diversity among stocked population 
from a same region and consequently caused an artificial 
increase of a regional based structuring. Breakdown of the 
genetic structure at various scales as also been observed in 
other studies where stocking strategies resulted in the 
homogenization of the genetic diversity among populations 
(Eldridge and Naish 2007; Perrier et al. 2011). 

Population and individual admixture following stocking 

Accordingly with the high genetic differentiation and the 
low proportion of shared diversity among unstocked lakes, 
most of the genetic clusters delineated by structure corre- 
sponded to a unique lake in the case of unstocked popula- 
tions. Almost no admixture was found in unstocked lakes 
in which almost all fish (97.0%) were assigned to the local 
clusters. These results reinforce the evidence that pro- 
nounced genetic drift and the very limited potential for 
migration among lakes have shaped the distribution of the 
contemporary genetic diversity in those lake trout popula- 
tions. These results are concordant with those reported by 
McCracken et al. (2013) in lake trout as well in other non- 
migratory salmonids (Koskinen et al. 2001; Ozerov et al. 
2010). 

While almost each unstocked lake corresponded to a 
unique genetic cluster with nearly no admixture, variable 
admixture rates were observed within stocked lakes as a 
function of stocking history. Indeed, most of the stocked 
lakes were highly admixed with their stocking source(s) 
and most of the individuals were assigned to nonlocal clus- 
ters or to the admixed category, while a smaller proportion 
of trout were assigned to the local cluster. These results 
demonstrate that stocking with nonlocal fish have caused a 
genetic displacement of the local gene pool. The occurrence 
of admixed individuals further indicates that stocked lake 
trout reproduced with local fish, potentially threatening the 
local genetic background through introgressive hybridiza- 
tion, as reported in other studies (Eldridge and Naish 2007; 
Miller et al. 2012). In other salmonids, admixture rates 
within stocked populations have also been found to be 
highly variable from almost absence of admixture up to a 
complete replacement of the native gene pool (Hansen and 
Mensberg 2009; Perrier et al. 2011). Interestingly, the 
genetic composition of stocked populations differed among 
administrative regions and can be related to the stocking 
modalities of each region. Indeed, in two of five regions, 
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stocking has been conducted less intensively with lower 
stocking intensity and a mix of local, admixed and nonlocal 
individuals were identified. In contrast, the three other 
regions were intensively stocked and stocking was con- 
ducted with high-stocking intensity and only nonlocal and 
admixed individuals were identified in stocked lakes. These 
results suggest that relatively large and intensive stocking 
may totally disrupt the original local genetic composition 
of wild stocked populations. 

We acknowledge that using an historical baseline prior 
to stocking may have allowed a better understanding of the 
temporal dynamic of admixture (Nielsen et al. 1999; Han- 
sen 2002; Finnegan and Stevens 2008; Metcalf et al. 2012; 
Perrier et al. 2013a,b), but such samples were not available 
for stocked lakes. Nevertheless, the genotyping of historical 
samples from some unstocked populations revealed high 
temporal stability over 40 years such as was found and dis- 
cussed for several salmonids (Tessier and Bernatchez 1999; 
Hansen et al. 2002; Vaha et al. 2008; Ozerov et al. 2010; 
Gow et al. 2011; Van Doornik et al. 2011) and thus rein- 
force our interpretation that admixture found in stocked 
populations most probably resulted from stocking rather 
than from naturally occurring changes in genetic composi- 
tion over time. Likewise, the high coherence between the 
well-documented stocking intensity and the observed 
admixture of stocked populations reinforces our conclu- 
sions. 

Stocking intensity largely explained admixture among 
stocked populations as populations' membership to the 
local cluster strongly decreased with the increase in stock- 
ing intensity. Stocked populations seemed to be highly sen- 
sitive to stocking with nonlocal fish as even small stocking 
intensities (eight fish/ha or a single stocking event) resulted 
in significant admixture and potential displacement of local 
gene pool. When stocking intensity increased up to 86 fish 
per hectare or up to 18 stocking events, all populations 
were subjected to a total displacement of the local gene 
pool and therefore to a complete and most probably irre- 
versible loss of the local genetic diversity. These results are 
in accordance with the few studies which documented 
stocking intensity to be correlated with the reduction of the 
genetic integrity in salmonids populations (Eldridge and 
Naish 2007; Hansen and Mensberg 2009; Marie et al. 2010; 
Perrier et al. 2013b). However, to our knowledge, this is 
the first study that clearly links admixture rates to stocking 
intensity on such a large geographic scale and a wide range 
of stocking modalities. Nevertheless, many factors may 
influence degrees of admixture between nonlocal and local 
fish such as the effective size of the recipient population, 
the environmental conditions and potential spatial or tem- 
poral reproductive isolation (Currat et al. 2008; Hansen 
and Mensberg 2009; Marie et al. 2012; Perrier et al. 2013b) 
and require further investigations. 



Admixture decreases following the cessation of stocking 
activities 

Few studies provided clear outcomes about the potential 
resilience of the original genetic diversity of wild popula- 
tions after the cessation of stocking (but see Perrier et al. 
2013b). Indeed, most studies compared genetic character- 
istics of the stocked populations between pre- and post- 
stocking samples, which generally highlighted an increase 
in admixture following stocking (Martinez et al. 2001; 
Hansen 2002; Susnik et al. 2004; Campos et al. 2007; 
Eldridge et al. 2009; Karaiskou et al. 2011; Pearse et al. 
2011; Glover et al. 2012; Perrier et al. 2013a). Here, while 
we could not compare samples collected prior or during 
stocking events, we found that the time elapsed since the 
last stocking event was a significant factor influencing the 
current admixture of populations. More precisely, we 
observed that after cessation of stocking, the genetic 
membership of stocked populations to their source popu- 
lations decreased over time to almost reach a null value. 
Different scenarios could be invoked to explain these 
results. One scenario could imply genetic drift, where 
ongoing divergence between source and stocked lakes as 
the last stocking event would increases with time. There- 
fore, the longer since the last stocking event, the higher 
would be the divergence between source and stocked 
populations. However, the strong genetic temporal stabil- 
ity we found for some unstocked populations is indica- 
tive of a low genetic drift in lake trout, which is in 
accordance with the species long generation time. An 
alternative scenario would assume that the observed 
decrease in genetic membership to source population 
reflects the purge of exogenous genes. This could be 
linked with a lower fitness of the stocked fish compared 
to the local ones (Araki et al. 2007a, 2008; Christie et al. 
2012; Milot et al. 2013). This could also imply that the 
impacts of stocking could be reversible if the exogenous 
genetic components are purged and that local individuals 
persist, as reported elsewhere (Nielsen et al. 1997; Hansen 
and Mensberg 2009; Perrier et al. 2013b). In the case of 
lake trout, such persistence of local individuals may be 
particularly long since lake trout is a long living species 
(up to 49 years). Potentially surviving old individuals 
may thus contribute to the conservation of the formal 
gene pool even in case of high-stocking intensity. The 
amount of time required to observe a complete resiliency 
could of course be variable depending notably on the 
stocking pressure, the relative amount of local fish per- 
sisting over time, potential introgressive hybridization, 
and fishing pressure (Evans and Willox 1991; Allendorf 
et al. 2001; Post 2013). However, in the absence of a his- 
torical baseline, we cannot confirm either a decrease in 
genetic membership to source populations, or an increase 
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in local genetic membership, as comparing historical and 
contemporary samples would have been necessary to test 
these hypotheses. We also acknowledge that these results 
could be obtained only for nine pairs of lakes with a 
low-stocking intensity of 17-22 fish per hectare. Evaluat- 
ing a wide range of different stocking intensities could 
have improved our understanding of resiliency processes, 
which could be different in heavily stocked lakes. 

Perspectives for management 

This study demonstrates that stocking profoundly altered 
the genetic integrity of wild lake trout populations, but 
that this alteration could be short-lived in some cases with 
relatively low stocking intensity. Our results on neutral 
genetic diversity showed genetic introgression to be 
strongly linked to stocking intensity, however, rates of 
introgression between different loci under selection could 
be highly variable as found for other stocked salmonids 
(Hansen et al. 2010; Meier et al. 2011; Lamaze et al. 2012, 
2013). Even though the extent of local adaptations in lake 
trout remains largely unknown, stocking could alter local 
adaptation in wild populations, as shown in other fish 
species (reviewed in Laikre et al. 2010). As preserving 
local adaptation should be of prime importance in conser- 
vation actions, it appears essential to further investigate 
the extent of local adaptation in lake trout and to quantify 
potential differences in introgression rates on genes which 
could confer adaptive advantage or disadvantage. Our 
results also emphasize the need to manage each lake trout 
populations independently and to adjust stocking strate- 
gies in order to conserve the genetic integrity and putative 
local adaptations of wild lake trout populations. Our 
results also provide fisheries' managers with new tools to 
define conservation targets for unstocked populations as 
well as for stocked populations. Indeed, a recovery strat- 
egy could be worth considering for those stocked popula- 
tions that did not reach a threshold value of total 
homogenization with their source, as recovery of the for- 
mal gene pool could be possible with time. In the case of 
the lake trout in the province of Quebec, results of this 
study have been used to propose conservation prioritiza- 
tion based on the genetic status of stocked populations. In 
addition with other criteria such as the habitat quality, 
fishing intensity and the population density, genetic rank- 
ing of populations has now been integrated in the deci- 
sion process to optimize stocking strategies in order to 
reduce the genetic impacts of stocking on wild lake trout 
populations. Conclusions of this study are clearly applica- 
ble for the species in all its distribution range as well as 
for other freshwater nonmigratory fish that might be 
highly structured, sensible to overexploitation and subject 
to stocking. 
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Additional Supporting Information may be found in the online version 
of this article: 

Table SI. List, characteristics and GenBank accession number of the 
markers used. 

Table S2. Pairwise F S j among populations, (spreadsheet). 
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Table S3. Pairwise F ST among the subset of lakes for which contem- 
porary and historical samples were available. Number in parentheses is 
the year samples were collected. 
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Table S4. If stocking occurred, the number of the source lake(s) are 
indicated in the 'Source(s) lake(s) used for stocking' column. A entry, 
indicates no stocking, (spreadsheet). 



644 



©201 4 The Authors. Evolutionary Applications published by John Wiley & Sons Ltd 7 (2014) 625-644 



